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ON THE STABLE RECOVERY OF A METRIC FROM 
THE HYPERBOLIC DN MAP WITH INCOMPLETE DATA 

PLAMEN STEPANOV, GUNTHER UHLMANN, AND ANDRAS VASY 


Abstract. We show that given two hyperbolic Dirichlet to Neumann maps associated to two 
Riemannian metrics of a Riemannian manifold with boundary which coincide near the boundary 
are close then the lens data of the two metrics is the same. As a consequence, we prove uniqueness 
of recovery a conformal factor (sound speed) locally under some conditions on the latter. 


1. Introduction 

Let {M,g) be a Riemannian manifold with smooth boundary. Let be the Laplace-Beltrami 
operator on M. In local coordinates, g is represented as a symmetric positive definite matrix 
g = {gij). The Laplace-Beltrami operator is given in local coordinates by 

= (det5)-|^(det5)^5*^^, 

where = {gij)~^, det^ = det{gij). Let A be a first order differential operator with smooth 
coefficients. Consider the following initial boundary value problem 

( (a| -Ag-A)u = 0 in (0, T) x M, 

(1) < u\t=o = dtu\t=o = 0 in M, 

[ '*^l(0,T)x9M = /) 

where T > 0 is fixed, / G R^^q^(( 0, T) x dM), and the subscript (0) indicates that the functions in 
that space vanish for t = 0. Denote by u = u{x) the unit outer normal to dM at x G dM. Define 
the hyperbolic Dirichlet-to-Neumann (DN) map Ag by 

Agf := 

where d^ is the exterior unit (in the metric g) normal derivative. The problem we study is if 
we can determine g, up to an isometry, in a stable way from the DN map. Uniqueness of such 
recovery has been established in [3], using the boundary control method developed by [2] that relies 
in the fundamental unique continuation result of Tataru EH. In this article we do not consider 
the determination of A. See [8] and the references there for results in this direction. Holder type 
stability estimates for a generic class of simple metrics were established in m- For the case of 
simple metrics with lower order terms stability estimates were given in [12] modulo the appropriate 
gauge transformations of the lower order terms, see also [Ij. We recall that simplicity of a metric 
means that the boundary is strictly convex and that the metric has no conjugate points. 

The goal of this note is to study stability for possible non-simple metrics assuming we know A 
locally. At present, we restrict our attention to recovery of the principal symbol which dictates 
the geometry. We show that in Theorem 3.1 that if the DN maps associated to to Riemannian 
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metrics are close enough in their natural norms, and g is known in a neighborhood of clM, then 
its lens relation is uniquely determined (exactly), and the problem is then reduced to uniqueness 
(not stability) of the lens rigidity one. We also prove a local version of this result. The lens 
rigidity problem, and the closely related boundary rigidity problem are well studied, see, e.g., 
[siEiiiiiiniiiaiiaiiaETiiso] and for some classes of metrics, g can be recovered up to isometry. 
The last two works deal with non-simple metrics. 

A result of a similar type for simple metrics was proved recently in [T]. In fact, they considered 
non-simple metrics under a certain condition for the conjugate points but it is not known if there 
is a metric for which this condition holds. In contrast to the earlier works however, it was noticed 
in [T] that closedness of the DN maps leads to exact equality of the metrics (up to an isometry). 
We explain the reason for this below. We want to point out that there are two factors which make 
such a result possible: the assumption that we know g near the boundary, and the fact that we 
work with the standard norm of the DN map. In fact, those results show that this norm may not 
be natural when studying stability questions. 

We also prove in section 4 new stability for metrics in the same conformal class using Theorem 
3.1 and the paper [20] of the authors. More precisely, if go is given, and g = c~‘^go, is in the same 
conformal class of go we want to recover c in a stable way. In this case there is no isometry involved. 
One of the reasons for considering the operator A is to include in the results the acoustic equation 
((9| — c^Agp) as well; then we can take g = c~^go but A is non-trivial. In this case c represents 
physically the sound speed of the medium. For simple metrics in the same conformal class Holder 
type stability estimates were proven in m- We remark is that in this case the result is for all the 
sound speeds not just a generic set. For other type of stability estimates for the sound speed see 
HI]. 


2. Preliminaries 

By PE], for any / G x dM), the problem ([T]) has unique solution u G C{[0,T] x 

n (^([O, r] X L‘^{M)) continuously depending on /; and moreover, 

(2) A : H'(iQ)((0,r) X dM) L^{{0,T) x dM) 

is continuous. Here, FI^^q^((0,T) x dM) is the closed subspace of H^{{0,T) x dM) consisting of 
functions vanishing at t = 0. One cannot expect that it would depend continuously on g and in 
particular on c when g = c~^go with go hxed. To illustrate this, consider the following well known 
simple example first: the shift operator Ucf{x) = f{x — c) with c a constant. This can be considered 
as the dynamics at time t = 1 associated to the transport operator dt + cdx- Then Uc is a unitary 
operator on T^(R) for any c but Uc is not a continuous family of operators in the operator norm. 
To see this, one constructs functions fn with norm one with support shrinking to 0. Then actually, 
\\Uci — Dcjll = 2 for Cl 7 ^ C 2 - In particular, ||C/ci — HcjH < 2 implies ci = C 2 - This can be extended 
to elliptic FIOs of order 1, for example, with a canonical relation a graph of a diffeomorphism. To 
prove a similar result, one constructs a sequence of distributions asymptotically concentrated at a 
single point in the phase space. In short, this is the idea of the proof of Theorem [1] below. 

Going back to A, it is an elliptic FIO of order one microlocally for covectors which are projections 
of characteristic covectors pointing strictly into R x M, and away from any neighborhood of t = 0, 
so we are in a similar situation. On the other hand, near t = 0 and x = y a dM, and directions 
pont strictly into M, A is a TDO of order 1 depending continuously on c. Near tangent directions 
however, the two Lagrangians intersect. The assumption that g is known near dM eliminates the 
need to study that intersection, where most of the difficulties lie. It also eliminates the need to 
study the diagonal (i.e., the transversal geodesics, at t = 0). The latter problem is well understood 
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in terms of uniqueness and stability, even for non-convex boundaries, see mm and the references 
there. 

We define the scattering relation next. Assume for convenience that dM is strictly convex, see 
also m for the more general case. Let d±SM be the open sets of all vectors (x, v) with x € 9M, 
V unit in the metric g, pointing outside/inside M. We define the lens relation 

L : d-SM d+SM 

in the following way: for each (x, v) G d^SM, L(x, v) = (y, w), where y is the exit point, and w the 
exit direction, if exist, of the maximal unit speed geodesic 'jx^v in the metric g, issued from (x,x). 
Let 

i : d-SM ^ RUoo 

be its length, possibly infinite. The metric g is called non-trapping, if i has a finite upper bound. 
In the latter case, we denote by Tq{M) the least one. 

It is convenient to parameterize d±SM by its projection on the unit ball bundle BdM of the 
boundary. In either case, the projection uniquely identifies one of the two unit vectors with the given 
projection by the choice of the sign ±. We will keep the notation L and i for the so parameterized 
quantities, as well. Given p' € BdM, the geodesic issued from p G d-SM having a projection p' 
will be also called the geodesic issued from p' and sometimes denoted by 7 p/. 

We will identify below vectors and covectors by the metric. 

3. Main results 

Theorem 1. Let fC (s d-SM be compact with maxyc-^ < oo. Let r± be two open subsets of dM so 
that 7r(K) C r_ andTToL{K) C r_|_. Let e > 0 andT — e > max^:-^. Then there exists 6 = S{JC) > 0 
so that if 

(3) 11^3 “ ^9llHi((0,£)xri)^.L2((o^'r)xr2) - 
for some other metric g, then Lg = Lg and ig = ig on K,. 

Proof. We will use solutions essentially supported near a single geodesic which may have conjugate 
points. In [T], this is done using Gaussian beams. We will use the semiclassical calculus with h > 0 
a small parameter, see e.g., [22]. The properties of the solutions then follow from the calculus, and 
there is no need to construct them explicitly even though we can: on each short enough interval 
for t, those solutions are given by the geometric optics construction with a suitable phase solving 
the eikonal equation, see e.g., [Hj. We still need to construct a reflected solution ui below but the 
latter needs to be constructed for a small time interval, and the same standard geometric optics 
construction then works well, see also m- 
Recall that the following functions 

F(z;zo,C°) = 

are called coherent states, for any fixed 7^ 0, see [22]. The semiclassical wave 

front set WF/i(T) consists of one point only, (zo,C°)- The extra factor h is chosen so that 1/C < 
||F||j:^i < C. Without that factor, F is unit in L^. 

Fix pq := (xq,^*^) G B*dM, and work in fixed local coordinates in some neighborhood Uq of xq. 
Set 

(4) /po = T(^,x;e/2,xo,-l,C°)• 

Here, z = (t, x). We can always assume that fp^ is localized near (e,xo) with a cut-off function. 
This is the coherent state on dM with WFhifp) at {to,xo,T^,f,^) = (eo/2,xo, —1,^*^). The standard 
geometric optics construction, see, e.g., [Ml [IS], implies that the solution u of O with / = fp^ 
has a semi-classical wave front on the zero bicharacteristic issued from where = 1 (the 
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norm is in the metric g), tf = in boundary normal coordinates, and is the positive 

root of other words, g^ is the outgoing (future pointing) unit covector with 

projection There are no singularities propagating back to t = 0 because of the zero Cauchy 
data there. If we replace tq = — 1 in ([1]) by r+ = +1, we get a similar singularity with r/)) < 0 
propagating along the bicharacteristic determined by — (in other words, by but for 

negative times). We are not going to use the second kind. 

Assume (as in the theorem) that the geodesic issued from po nt t = e/2 is non-trapping, i.e., 
hits dM again in r 2 , transversely, for a finite time ti = i{po), at (ti,xi). Then the wave front set 
of u is on the bicharacteristic over that geodesic (see, e.g., [22]), until it hits dM, again, then it 
reflects (see m, etc. Let for a moment uq be as u but propagating outside of M near (ti, xi); with 
some smooth extension of g there. The map / i-)- 5 i/Uo|rx 9M) for any / with wave front set near 
po (notice that tangential directions are avoided) is an /i-FIO of order 1 with a canonical relation a 
diffeomorphism. Let ui be the reflected solution, i.e., uq + rti = 0 on dM near {to,xo), and ui has 
wave front set obtained by that of uq by reflection, see usms]. Then in M, in some neighborhood 
Ui of (ti, xi) in R X dM, u is given hy u = uq + ui and the map / i-)- d^u\u 2 is an /i-FIO of order 1 
again, with principal symbol equal to twice that of the FIO above. Therefore, if we choose Ti > ti 
so that it is less than the time t 2 for which the second reflection takes place, WF^ ('U[o,Ti]xPA 4 ) lies 
over UoUUi. 

So far po = (xo,p^) was fixed and the construction is not invariant under change of variables. 
Let us now vary p over the whole 1C. For p near each such po, we have an uniform estimate from 
below for ||A/p||^ 2 ([;^) for 0 < /i < /iq with some /iq > 0 with Ui = Ui{xo,g^) as above. Then we 
can chose a finite sub-cover with the same property and define fp globally by a partition of unity. 
Therefore, there exists (5o > 0 so that 

(5) 2(5o < ||A/p||j^2((£ 'r)x9M)) Vp €/C, 0 < h < ho. 

Notice that above, we excluded the interval (0,e). 

Let u be constructed as above, but related to g. Under the assumptions of the theorem, using 
finite speed of propagation, u{t,x) = u{t,x) for t near e/2 and all x. Then the normal derivatives 
near (e/2,xo) equal as well, i.e., A^/p = A^/p there. Assume that {Lg,ig) and {Lg,£g) do not 
coincide in 1C, which includes the possibility that for some p G K., 7 p, related to g, is trapping, 
i.e., ig{p) = oo. Then there exists p G U so that {Lg{p),ig{p)) ^ {Lg{p),ip{p)), or the latter is 
undefined. Then near WF/j(Ap/p), we have Agfp = 0{h°°) and over (e,T) x dM, A.gfp and Agfp 
have non-intersecting wave front sets if T is chosen larger than e/2 -|- i{p) but smaller than the 
second time the reflected ray hits dM again. Then they are orthogonal modulo 0{h°°). Then 

, 11(^9 “ ^9)/plli2((ro,T)xaM) = ll^9/plli2((£,T)xPM) + ll^9/plli2((£,T)X(9M) + 0{h°°) 

^ Il^<7/plli2((£,T)xPM) -0{h°°) 

This however contradicts m and dSj) for 0 < /i < ho{p,g) -Cl. □ 

In particular, with Fi = F 2 = dM, and if p — p is a priori supported away from dM, then 5 <C 1 
in (|3|) implies Lg = Lg and ig = ig globally, which is one of the results in [T] . 

Note that the l.h.s. of (|3|) in Theorem [1] can be microlocalized by replacing F± with open sets 
in the phase space BdM] then we remove the projection tt in the formulation of Theorem [H To 
localize there, we can choose suitable zero order TDOs with principal symbols supported in those 
sets, and use them as cutoffs. 


4. Recovery of a sound speed 

To get results for non-simple metrics and partial data, we apply a recent result obtained by the 
authors in [20]. Assume in this section, that p = c“^po with po fixed. We showed in [20], that one 
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can recover c near some p € dM in N, if g is given, by the lens relation known locally near SpdM. 
Then we extended this local result globally, assuming the foliation condition below. Note that we 
can have conjugate points under the foliation condition, without a restriction of their type. Also, 
the global arguments in [20] are based on layer stripping arguments, and work even if recovery of 
the wave front set from the X-ray transform may not be possible locally. Below, M is an open 
extension of M. 

Definition 1. Let Mq G M be compact. We say that Mq can he foliated by strictly convex hyper¬ 
surfaces if there exists a smooth funetion p : M —>■ [0, oo) which level sets St = p~^{s), 0 < s < S 
with some S > 0, restricted to M, are strictly convex viewed from p“^(0, s) for g; dp / 0 on these 
level sets, Sq n M = 0, and Mq C p~^(0, S]. 

Examples of foliations are the geodesic balls for metrics of negative curvature, or the surfaces 
{x G M\ dist(x, dM) = e}, for e -C 1 if dM is strictly convex. Any Mq as above is non-trapping in 
the sense that any maximal geodesic in Mq is of finite length. We denote the least upper bound of 
those geodesics by Tq{Mq). 

Then we get the following local result, as consequence of Theorem 5.2 in [20]. Note that we use 
the uniqueness part of the latter only. 

Corollary 1. Let n = dimM > 3. Assume that Mq C M can he foliated by strictly convex 
hypersurfaces p~^{s), 0 < s < S'. Let T C dM he a neighborhood o/Tq := p“^(0,S) AdM. Let 
c = c in some neighborhood o/Tq in M. Then if ^ holds with g = c~‘^gQ, g = c~‘^gQ, and with 
Ti = r 2 = T, and Tq > Tq{M), we have 

(7) c = c on Mq. 

Clearly, we can require above only the geodesics hitting Tq to be non-trapped, with Tq being an 
upper limit of their lengths. 

Combining Theorem [1] with [20], we immediately get the following. 

Corollary 2. Let n = dimM > 3. Assume that M = Mq U Mi is relatively open and can be 
foliated with strictly convex surfaces with respect to the metric g, and Mi is simple. Assume, as 
above, that c = c near dM. If ^ holds with Ti = r 2 = dM, T > Tq{M), then c = c. 
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